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On the Theory of the Capillary Tube, 
By Lord Eayleigh, O.M., F.li.S. 

(Received October 26, 1915). 

A recent paper by Richards and Coombs* discusses in some detail the 
determination of surface tension by the rise of the liquid in capillary tubes, 
and reflects mildly upon the inadequate assistance afforded by mathematics. 
It is true that no complete analytical solution of the problem can be obtained, 
€ven when the tube is accurately cylindrical. We may have recourse 
to graphical constructions, or to numerical calculations by the method of 
Runge,f who took an example from this very problem. But for experimental 
purposes all that is really needed is a sufficiently approximate treatment of 
the two extreme cases of a narrow and of a wide tube. The former question 
was successfully attacked by Poisson, whose final formula [ (18) below] would 
meet all ordinary requirements. Unfortunately doubts have been thrown 
upon the correctness of Poisson's results, especially by Mathieu,$ who rejects 
them altogether in the only case of much importance, i.e. when the liquid 
wets the walls of the tube — a matter which will be further considered later 
on. Mathieu also reproaches Poisson's investigation as implying two different 
values of h, of which the second is really only an improvement upon the 
first, arising from a further approximation. It must be admitted, however, 
that the problem is a delicate one, and that Poisson's explanation at a critical 
point leaves something to be desired. In the investigation which follows I 
hope to have succeeded in carrying the approximation a stage beyond that 
reached by Poisson. 

In the theory of narrow tubes the lower level from which the height of 
the meniscus is reckoned is the free plane level. In experiment, the lower 
level is usually that of the liquid in a wide tube connected below with the 
narrow one, and the question arises how wide this tube needs to be in order 
that the inner part of the meniscus may be nearly enough plane. Careful 
•experiments by Richards and Coombs led to the conclusion that in the case 
•of water the diameter of the wide tube should exceed 33 mm., and that 
probably 38 mm. suffices. Sucli smaller diameters as are often employed 
(20 mm.) involve very appreciable error. Here, again, we should naturally 
look to mathematics to supply the desired information. The case of a straight 

* ' Journ. Amer. Cbem. Soc.,' No. 7, July, 1915. 

+ ' Math. Ann.,' vol. 46, p. 175 (1895). 

J ' Theorie de la Capillarit6,' Paris, 1883, pp. 46-49, 
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wall, making the problem two-dimensional, is easy,* but that of the circular 
wall is much more complicated. 

Some drawings (from theory) given by Kelvin, figs. 24, 26, 28,t indicate 
clearly that diameters of 1*8 cm. and 2*6 cm. are quite inadequate. I have 
attempted below an analytical solution, based upon the assumption that the 
necessary diameter is large, as it will be, if the prescribed error at the axis 
is small enough. Although this assumption is scarcely justified in practice, 
the calculation indicates that a diameter of 4*7 cm. may not be too large. 

As Eichards and Coombs remark, the observed curvature of the lower part 
of the meniscus may be used as a test. Theory shows that there should be 
no sensible departure from straightness over a length of about 1 cm. 

The Narroiu Tube. 

For the surface of liquid standing in a vertical tube of circular section, we 

have 

X dzjdx IP . 

in which z is the vertical co-ordinate measured upwards from the free plane 
level, X is the horizontal co-ordinate measured from the axis, i/r is the angle 
the tangent at any point makes with the horizontal, and ' a^ = Tgp,X where 
T is the surface-tension, g the acceleration of gravity, and p the density of 
the fluid. The equation expresses the equilibrium of the cylinder of liquid 
of radius x. At the wall, where ^ = r, yfr assumes a given value (^Tr — i), 
and (1) becomes 

ah" cos i =: \ zx dx. (2) 

Jo 

If the radius (r) of the tube is small, the total curvature is nearly 
constant, that is, the surface is nearly spherical. We take 

zz=l^^(c^-.x^)-]-u, (3) 

where / is the height of the centre and c the radius of the sphere, while u 
represents the correction required for a closer approximation. If we omit u 
altogether, (2) gives 

■^ Compare ' Phil. Mag,,' vol. 34, p. 309, Appendix, 1892 ; ' Scientific Papers,' vol. 4, 
p. 13. 

t The reference is given below. 

I It may be remarked that a^ is sometimes taken to denote the double of the above 
quantity. 
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Also, if h be the height at the lowest point of the nienisciis, the quantity 

directly measured in experiment, 

h^l~c, (5) 

In this approximation t\c — cos i, and thus in terms of o 

ah'^jc = \r\h-{-G)-{-^^{c^ — T^Y!^'—^c^. (6) 

When the angle of contact {i) is zero, c = r, and 

a2 = ir(/^ + ir), (7) 

the well known formula. 

When we include ii, it becomes a question whether we should retain the 

value of c, i.e. r sec i, appropriate when the surface is supposed to be exactly 

spherical. It appears, however, to be desirable, if not necessary, to leave the 

precise value of c open. Substituting the value of z from (3) in (1), we get, 

with neglect of {d/ujdxf, 



1 + 



du {g^— x^y^^ 3( c^ — a)^y Idu \^ ___ c Vlx^ (c^ — x^^f^ — c^ 

clx xc^' 2c* ^dx cC^x^ 



T + - 3 + 



'X 



ivX LvX 

Jo 



• (8) 



For the purposes of the next approximation we may omit {dujdxy and the 
integral, which is to be divided by ci?. Thus 

dih __l cl ___-|\ c^x ^ c^ 1 , ^^ /Q\ 

dx \2a^ /(c^— ^^)^/^ oa? x{c^—x^fl'^ "da^x' 

and on integration 



It 



'to 

cl C^ ^ \ 1 . C^ 



^iirb-V^TW^^i^''^^'''^^^-''^^-'''- ^''^ 



We suppose with Poisson and Mathieu that 



'^3-A.-^ = ^> (11) 



2a^ 3a^ 



c^ 



so that '^ = -5-2 log {c + ^(c^— a?^)} + 0, (12) 



Set 

mg to = ^ ^ ^. (13) 

dx da^^ x^/{&^—x^) 

To determine c we have the boundary condition 

, . _ fdz \ _ T , /du\ 

- /l-^ Z-v/Cf-^") ! (14) 



which gives c in terms of i and r. Explicitly 



r 7'^ sin^ i 



cos t 00/^ cos^ t ( 1 4- sm %) 



These latter equations are given by Mathieu. 
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We have now to iBnd the value of a^ to the corresponding approximation. 
Eor the observed height of the meniscus 



^3 

h = l—c + U:,=.o = l—c + G'i-^-^log(2c); (16) 

OC{/ 

and 

ah^ cos i = zxclx =^ — {l-{'G)'\--^ Uc^ + r^^l^ — c^} + {u— G)x dx 
Jo ^ ^ 



«/ 







= (^i + c)|' + |{(c^-^T'-«'} 



C3 



3a^ 



"'^ log£+4£!=:!3 + i^ + v'(f-^'')}'-V(g^-r^)1. (17) 



In the important case where i = 0, the liquid ivettmg the walls of the 
.tube, c = r simply, and 

«"='2(*+i)-6T'("«^-5: 

= -|r(A+ ^/•-01288r7A). (18) 

This is the formula given long since by Poisson,* the only difference being 
;that his a^ is the double of the quantity here so denoted. 

It is remarkable that Mathieu rejects the above equations as applicable to 
the case i = 0, e = r, on the ground that then dujdx in (13) becomes infinite 
when X = r. But d^{T'^—x^)ldx, with which dujdx comes into comparison, 
:is infinite at the same time ; and, in fact, both 

-—■{t^—x^Yi^ and -— • (7-^—o(?f, 

dx \dx/ 

in equation (8) vanish when x = r. It is this circumstance which really 
determines the choice of /in (11). 

We may now proceed to a yet closer approximation, introducing approxi- 
mate values of the terms previously ueglected altogether. From (13) 

^ ^ \dxl 9a^x^ ^ ^ ^ ^ ^ ^ ^ ^' 

.and from (12) 

toxdx = ^Gx^+ w-^ [x^ log {c+^(c^— ^^)} + 2^*^'~^\/(^^'~^^) + iC^"""^)]' 



■^ 



' Nouvelle Theorie de 1' Action Capillaire,' 1831, p. 112. 
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Thus 



du _ /c(Z + C) 



dx 



2w 



2 



1 



C^X 









(c^ — x^)^l'^ 3 c(?x (c^— x^y^'^ 3 a^x 



%xa^{c^—x^Y^^ 



Za 



4.(c2_a;2)2_2c(c2_a;2)3/3 



+ A;2 log { c + v/ (c^ - »'") } + ij - c V(c' - *3') 



(19) 



'Z6 = 



c([+C)__ c 



O 9 

36:r 



1 



"1 c^ __ 



Fx + ^-^^g{^ + \/(^^-^')} 



c 



5 "- 



6a* L 



-2 log {c + v'Cg'— ^')} +'^^"'""'^'^-1 



c 



2. /r^2_^2) 



$v(^^ 



<y/((3^ — 0:^^) 



log{c + ^(G2~a;2)} 



constant. 



We have now to choose L or rather (Z + C), and it may appear at first sight 
as though we might take it almost at pleasure. But this is not the case, at 
any rate if we wish our results to he applieahle when c — r. For this purpose 
it is necessary that {duldx)rX{T'^—x^) he a small quantity, and only a 
])articular choice of (/ + C) will make it so. For w^hen x ~ c ~ r, 



dxjr 



■X" 



'r{l-\-G) 



1-1 <a 



^y{T^'-X^) L 2a 



2 



1 __- _!__. -|- 



3a2 ' 6an ^ ^2 



IM 



%o} 



+ terms vanishing when x = r. 
We must therefore take 



c([+C) 

la' 



v2 



C^ I 



?)a^ 6(2^ 



1\ 



i"l;2+^(^og^ + o) = ^^ 



(20) 



making 




lb = 


C3 / C2 




6a^ 



^^ho- ^"^^^^^^-"'^^^ 



G 









+ G\ 



(21 



It should he noticed that ii so determined does not become infinite when 

c — r and x == r. For we have 



ivv — — 



^' 



■^O 



rA 



+ 



Ga'^ t)a'^.\/(7^^ — ^' 



-l0g|l+V^^'''^ ^'H+C^ r=. C. 



Also with the general value of c 



Uo 



G 



.3 



3a 



2 



1 



.2 



261^ 



., lo^ 2 + C^ 



(22) 



As before 



h = l—c-i-tiQ, 
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and 



rrr cos 



. ^ (i±^.,(lz:ff=l^ [\u-0').c^ 



m 



Jo 



ta 



^■^'-uA^-2^r-'^^-^ 



(c^ ~~r^yi^ — ^3 



fr 






{U'^G')xdx. (23) 







The integral in (23) can be expressed. 
We find 



i 






L2 



c 



.5 r 






6a* 



VS ,..__ //'S . n',2'^12 ma 



3c 



2 



i{c + ^(c2-r2)} /logl±.i/^^!i:^^-lT 



+ 2c2(log2-l) 



(24) 



The expression for ^'a^cosi in terms of c is complicated, and so is the 
relation between c and i demanded by the boundary condition 

'd'ii' 



cott 



r 



"X" 



^(C^ — r^) \dx jx^r' 



(25) 



But in the particular case of greatest interest {% = 0) much simplification 
ensues. It follows easily from (25) that c = r. When we introduce this 
condition into (24), we get 



{v, — G')xdxz= T7^+s7i(21og.2— ij). 







12*2 ' 6a* 



(26) 



and accordingly 



a-=|(A+0 



5r^ 



^/log24) + 3^,(31og2^2).. 



(27) 



Hence by successive approximations 



a 



T 

2 



3 



'-^^^+-^-^(log2-i)+j^^(321og2-21) \ 



r 



,3 



} 



(28) 



2r2 

1 iOS" 25 — 1 + 

= |r{A + ir^01288r7>^ + 01312rVA2}, 

If the ratio of r to ^ is at all such as should be employed in experiment, this 
formula will yield a^, viz., T/^p, with abundant accuracy. 

Our equations give for the whole height of the meniscus in the case ^ = 0, 
e = r, 

(29) 



,_,„ = ,^.,,_,„ = ,__^l__^log2. 

Another method of calculating the correction for a small tube, originating 
apparently with Hagen and Desains, is to assume an elliptical form of surface 
in place of the circular, the minor axis of the ellipse being vertical. In any 

VOL. xcii. — A. Q 



■r 



o hr , hr / L , 4:r\ r f . , r r- , lt \ 
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cas e this should allow of a closer approximation, and drawings made for 
Kelvin* by Prof. Perry suggest that the representation is really a good one. 

If the semi-axis minor of the ellipse be yS, the curvature at the end of this 
a xis is /3/r^, and in our previous notation /3 = hr'^j^a?. Also, i being equal 
to 0, 

and a? = ir(/i + -|-/3) = \}it{\^-t''\U^). (30) 

Th is yields a quadratic in a^ ; hence 

3 yA m^j 

= ir{7i + -i-r-0-llll rVA + 0-0741 rVA^} (31) 

approximately. It will be seen that this differs but little numerically from 
(28), which, however, professes to be the accurate result so far as the term in 
r'^jh? inclusive. 

The Wide Tube, 

The equation of the second order for the surface of the liquid, assumed to 
be of revolution about the axis of z, is well known and may be derived from 
(1) by differentiation. It is 

dx^ xdx\ \dx) J a?\ \dxj J 

If dzjdx be small, (32) becomes approximately 

<i% 1^_ z __ ?>z (dzY_l(dzY /oQ\ 

dx^ xdx a^ 2a^\dxj x\dxj' 

In the interior part of the surface under consideration {dzjdxY may be 
neglected, and the approximate solution is 



(AJ , iAJ 



z = /ioJo(W^O = /^olo(^/f0 = ^^0 j^l + 2^ + ^2 42 a^ "^"j ^ ^^^^ 

Jo denoting, as usual, the BesseFs, or rather Fourier's, function of zero order, 
and 7^0 being the elevation at the axis above the free absolutely plane level. 
For the present purpose li^\^ to be so small as to be negligible in experiment, 
and the question is how large must r be. 

When 7^0 is small enough, x/a may be large while dzjdx still remains small. 
Eventually dzjdx increases so that the formula fails. But when ^ is, large 
enough before this occurs, we may if necessary carry on with the two- 
dimensional solution properly adjusted to fit, as will be further explained 

^ ' Proc. Roy. Inst.,' 1886 ; ' Popular Lectures and Addresses,' I, p. 40. 
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later. In the meantime it will be convenient to give some numerical examples 
of the increase in dzjdx. In the usual notation 

'^=hiJA (35) 

ax a \al 

and the values of Ii, up to xja = 6, are tabulated.* 

In the case of water a = 0*27 cm. If we take IiQfa = O'Ol, and x/a = 4^ 
we have dz/dx = 0*098, so that {dzjdxf is still fairly small. Here for water 
liQ = 0'0027 cm. and 2x — 2*2 cm. A diameter of 2*2 cm. is thus quite 
insufficient, unless an error exceeding 0*003 cm. be admissible, i^gain^ 
suppose haja = 0*001, and take x/a = 6. Then dz/dx = 0*061, again smalL 
For water ho — 0*00027 cm., and 2x ■==. 3*2 cm. This last value of \ is about 
that (0*003 mm.) given by Eichards and Coombs as the maximum admissible 
error of reading, and we may conclude that a diameter of 3*2 cm. is quite 
inadequate to take advantage of this degree of refinement. 

We may go further in this example without too great a loss of accuracy. 
Eetaining h^/a = 0*001, let us make x/a =7. I find Ii (7) = 156, about, 
so that the extreme value of dz/dx is 0*156, still moderately small Here 
2x = 3*8 cm., which is thus shown to be inadequate in the case of water. 

But apart from the question of the necessary diameter of tube, information 
sufficient for experimental purposes can be derived in another manner. The 
initial value of z (on the axis) is ho ; and z = 2ho when Io(x/a) == 2, ie. when 
X = I'Sa. FoT the best work Jiq should be on the limit of what can be 
detected and then ho and 2ho could just be distinguished. The observer may 
be satisfied if no difference of level can be seen over the range x = + 1*8^ ;. 
in the case of water this range is 2 x 1*8 x 0*27 = 0*97 cm., or say 1 cm. 

It has already been remarked that when ho is small enough x/a may 
become great within the limits of application of (35). To shorten our 
expressions we will take a temporarily as the unit of length. Then when x 
is very great. 



e^ 



^^^"> = ^»^"> = 7(2^- (^^> 



^=tant=^ = -7^^; (37) 



Thus if i/r be the angle the tangent to the curve makes with the horizontal, 

dz _ hoe^ 
dx ^{27rxy 

an equation which may be employed when ho is so small that a large x is 
consistent with a small i/r. 

In order to follow the curve further, up to i/r = -|7r, we may employ 
the two-dimensional solution, the assumption being that the region of 

■^ Brit. Assoc. Kep. for 1889 ' ; or Gray and Mathews' ' Bessel's Functions,' Table VI. 

Q 2 
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moderate ^ occupies a range of x small in comparison with its actual value, 
i.e. a value not much less than r, the radius of the tube. On account of the 
magnitude of x we have only the one curvature to deal with. For this 
curvature 

l=^ = ^sm^ = ^, (38) 

XV Wo Ujty 

SO that \'^ = C— cos-^lr =: 1— cos'\|r, 

since when i|r = 0, ;2;^ is exceedingly small. Accordingly 

= 2 sin l'^. (39) 

Also dx = - — ?-™ == [- — -- — 2 sin I -iir ) dik-^X 

tan'»|r xsml-^lr ^ ^ j ^^^ ^ 

and X = log tan (J-^) + 2 cob \'^ + G^ (40) 

The constant is determined by the consideration that at the wall (x = r), 
^ =z= I TT ; thus 

T—x = log tan (7r/8) + '\/2--log tan (^ilr)— 2 cos (|i|r) 

= log tan (tt/S) + v^2 - 2 + 2 log 2 - log -^^ (41) 

since '^ is small. 

The value of x is supposed to be the same here as in (37), so that 

X— log -^ + 1 log (27ra?) — log /i/o, (42) 

whence on elimination of '^ and restoration of a, 

TJa = -log (v/2 + 1)4-^^2-2 + 2 log 2 + -I log (27ra3/a) -log (^oM (43) 

With sufficient approximation, when Ao is small enough, we may here 
substitute r for x, and thus 

r/a—^log(r/a) = — log(v/2 + l) + \/2 — 2 + 21og2 + |log(27r)— log(feo/<^) 

= 0-8381+ log (a//io). (44) 

This formula should give the relation between r/a and h/a when h/a is 
small enough, but it is only roughly applicable to the case of greatest 
interest, where a/7io = 1000, corresponding to the accuracy of reading found 
by Eichards and Coombs. In this case 

0-8381 + log (a//io)= 7-H6. 

For this value of r/a we should have |log(r/«^) = 1-024. It is true that 
according to (44) r/a will be somewhat greater, but on the other hand the 
proper value of x (replaced by r) is less than r. We may fairly take 

r/a = 7-746 + 1-024= 8-770, 

making with a = 0*27 cm, 

2r = 4-74 cm. 

This calculation indicates that a diameter greater even than those con- 
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templated by Eichards and Coombs may be necessary to reduce Iiq to 
negligibility , but it must be admitted that it is too rough to inspire great 
confidence in the close accuracy of the final number. Probably it would be 
feasible to continue the approximation, employing an approximate value for 
the second curvature in place of neglecting it altogether. But although the 
integration can be effected, the work is rather long. 

l^Aclded November 17. — Since this paper was communicated, I have been 
surprised to find that the problem of the last paragraphs was treated long 
ago by Laplace in the ' Mecanique Celeste '* by a similar method, and with 
a result equivalent to that (44) arrived at above for the relation between 
the radius of a wide tube and the small elevation at the axis. Laplace uses 
the definite integral expression for lo, and obtains the approximate form 
appropriate to large arguments. In view of Laplace's result, I have been 
tempted to carry the approximation further, as suggested already. 

In the previous notation, the differential equation of the surface may be 
written 

?HL±i±4.E5i± = ^. (45) 

dz X a?' 

In the first approximation, where the second curvature on the left is 
omitted, we get 



2a^ 



1— cos T|r = 2 sin^-^, 



zq being the elevation at the axis, where i|r = 0. For the present purpose 

zq^ is to be regarded as exceedingly small, so that we may take at this stage, 

as in (39), 

z=^2a sin (|-^). (46) 

We now introduce an approximate value for the second curvature in (45), 
writing ^ = r, where r is the radius of the tube, and making, according 
to (46), 



z / / w z^ 



^^ = W^-i^- ^''^ 



On integration 



C^cos^ = ^, + ^(l^f^'^' = il + |^cos4, (48) 

^ 2a^ '6r\ 4aV 2a^ 3r 2 ^ ^ 

on substitution in the small term of the approximate value of z. When 
'\jr = 0, z^ is very small, so that C = 4:a/S7% and 

Z o • '^ . 2al— COS^(i'\|r) ,.^v 

- = 2 sm -x + . ,^y\ (49) 

a 2 or smji/r 

is the second approximation to z. 

^ Supplement au X® Livre, pp. 60-64, 1805, 
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From (49) 

We are now in a position to find x by the relation 

X = ^COt^Jr{dz/d'\fr)dyJr, (51) 

the constant of integration being determined by the correspondence of ^ = r, 
i|r = I TT. Thus 

^ = log(v/2---l) + ^2 + ^|--2 + ^^+|log(v/2-l)j--logtan| 

-^2cosf-fii|-^ ^ ^ - + 2Bin4^glogl r:^^f^^ T, (62) 
2 3r L2(l~cos-|f) 2 2^ sinj-^ J ^ ^ 

giving when 'yjr is small 

IZlE =, a^^ -(l + ^logir, (53) 

a 3r \ 2r/ ^ ^ ^ 

where ^ = log(v/2-l) + v/2 + log4-2 = - 0*0809, (54) 

—1-^ = log2 + Jlog(v/'2-~l) + iv'2-7/12 = -0-0952. (55) 

The other equation, derived from the flat part of the surface, is 

dx a \a/ ^(27Tx/a) \ SxJ' 
in which x/a is regarded as large ; or 

E=iog^ + log£+-|log?^+|. (57) 

In equations (53), (57) x and ^lr are to be identified. On elimination of -^ 
r -, €c oL — aS/Sr , r—x , -. , 2'nx , Sa ,p;o^ 

«-^'^'/v=T+5t + 2^ + ^^^°^— + 85' (^^^ 

in which we may put 

log = log hlog 1 = log — 



Thus 



a ■ a \ T I a r 

3a _ 3ce/^ I '^~"^\ 
8a? 8r \ r J ' 

— logy- = ■ , , Vo +-Vt-^+ -5" +i ^-Qg — ' (^^) 

in which, since ^ is nearly equal to r, a(r — x)/ST^ may usually be neglected. 
Also, in view of the smallness of a and /3, it is scarcely necessary to retain 
the denominator l-{-a/27% so that we may write 

5^0 



— logf-= -0-0809 + 0-2798~ + ilok — 



0-8381 + 0-2798 a/r + J log {tjo). (60) 
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(61) 



(62) 



The effect of the second approximation is the introduction of the second 
term on the right of (60). 

To take an example, let us suppose as before that a/ho = 1000, so that 
log (a/h(^) = 6*908. By successive approximation we find from (60) 

r/a = 8*689, 

so that if a = 0*27 cm. (as for water), 

2r = 4*79 cm. 

The correction to Laplace's formula is here unimportant. 

The above is the diameter of tube required to render ho negligible 
according to the standard adopted. 

It may sometimes be convenient to invert the calculation, and deduce the 
value of ho from the diameter of the tube (not much less than 4 cm.) and 
an approximate value of a. For this purpose we may use (60), or preferably 
(59), taking x = fr for instance. The calculated value of h^ would then be 
used as a correction. The accompanying small Table may be useful for this 
purpose. 



r/a. 


~logio(V»)- 


Difference. 


hQJa, 


6 

7 

8 

9 

10 


1 -8275 

2 -2319 

2 -6399 

3 -0508 
3 -4639 


-4044 
-4080 
-4109 
-4131 


-0149 
-0059 
0-0023 
-00089 
-00034 



We have supposed throughout that the liquid surface is symmetrical about 
the axis, as happens when the section of the containing tube is circular. It 
may be worth remarking that without any restriction to symmetry the 
differential equation of the nearly flat parts of a large surface may be taken 
to be 



so that z may be expressed by the series 



z 



a 



2 



0, 



(63) 



z = Ao lo{r/a) + (Ai cos ^ + Bi sin 0) Ii (r/a) 

+ (A3 cos 20 + B2 sin 20) hir/a) + ... 
r, denoting the usual polar co-ordinates in the horizontal plane.] 



(64) 



